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A SIGNAL PROCESSING METHOD FOR DETERMINATION OF A PARAMETER 
OF A SYSTEM GENERATING THE SIGNAL 

The present invention relates to a method for determination 
5 of a parameter of a system generating a signal containing 
information about the parameter. 

The method may be used for identification of sound or speech 
signals, such as in speech recognition, or for quality 

10 measurement of audio products or systems, such as 

loudspeakers, hearing aids, telecommunication systems, or for 
quality measurement of acoustic conditions. The method of the 
present invention may also be used in connection with speech 
compression and decompression in narrow band 

15 telecommunication . 

The method may also be used in analysis of mechanical 
vibrations generated by a manufactured device during ope- 
ration e.g. for detection of malfunction of the device. 

20 

The method may further be used in electrobiology for example 
for analysis of neuroelectrical signals such as analysis of 
signals from an electroencephalograph, an electromyograph, 
etc . 

25 

BACKGROUND OF THE INVENTION 

Prior art methods of signal processing are based on a short 
time Fourier transform of signals and it is assumed that the 
30 signals are steady state signals. 

In steady state analysis the signal is assumed stable in the 
period the signal is analysed and the steady state spectrum 
is calculated. 

35 

In real life steady state signals do not occur and steady 
state analysis does not provide sufficient knowledge of 
phenomena within various scientific and technological fields. 
Consider for example speech analysis. The human ear has the 
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ability to simultaneously catch fast sound signals, detect 
sound frequencies with great accuracy and differentiate 
between sound signals in complicated sound environments. For 
instance it is possible to understand what a singer is 
5 singing in an accompaniment of musical instruments. 

It is assumed that the cochlea in the human ear can be 
regarded as an infinite number of band-pass filters within 
the frequency range of the human ear. 

10 

The time response f (t) for one band-pass filter due to an 
excitation can be separated into two components, the 
transient response, fc{t), and the steady state response, 

f3(t), f (t)=ft {t)+f3(t) , 

15 

Traditional signal processing is based on the steady state 
response fs(t), and the transient response ft (t) is assumed 
to vanish very fast and to be without importance for the 
perception, see for example "Principles of Circuit 
20 Synthesis", McGraw-Hill 1959, Ernest 5. Kuh and Donald O. 
Pederson, page 12, lines 9-15, where it is stated that: 

"only the forced response is considered while the response 
due to the initial state of the network is ignored" . 

25 

Thus, when students are introduced to the world of signal 
analysis, they learn that the transient response, i.e. the 
response due to the initial state of the network should be 
ignored because it vanishes within a very short period of 
30 time. Furthermore, it is rather difficult to analyse these 
transient signals by use of traditional linear methods of 
analysis. 

The ability of the human ear to hear very short sounds and at 
35 the same time detect frequencies with great accuracy is in 
conflict with the traditional filterbased spectrum analysis. 
The time window (twice the rise time) of a band-pass filter 
is inversely proportional to the bandwidth, tw=2/ ( f u-f i ) , 
where fi is the lower cut-off frequency and fi, is the upper 
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cut-off frequency. 

ThuS/ if a rise time of 5 ms is required the consequence is 
that the frequency resolution is no better than 400 Hz. 

5 

As the detection of these transients is in conflict with a 
high frequency resolution, the detecting by the himan ear of 
these transients must take place in an alternative manner. It 
has not been examined how the human ear is able to detect 

10 these signals, but it might be possible that the cochlea, 

when no sounds are received, is in a position of rest, where 
the cochlea will be very broad-banded. When a sound signal is 
received, the cochlea may start to lock itself to the 
frequency component or components within the signal. Thus, 

15 the cochlea may be broad-banded in its starting position, but 
if one or more stable frequencies are received the cochlea 
may lock itself to this frequency or these frequencies with a 
high accuracy. 

20 Today it is known that the nerve pulses launched from the 
cochlea are synchronized to the frequency of a tone if the 
frequency is less than about 1.4 kHz. If the frequency is 
higher than 1.4 kHz the pulses are launched randomly and less 
than once per cycle of the frequency. 

25 

Signal processing based on filter bank spectrum analysis is 
disclosed in GB 2213623 which describes a system for phoneme 
recognition. This system comprises detecting means for 
detecting transient parts of a voice signal, where the 

30 principal object of the transient detection is the detection 
of a point where the speech spectrum varies most sharply, 
namely, a peak point. The detection of the peak points is 
used for a more precise phoneme segmentation. The transient 
analysis of GB 2213623 is based on a spectrum analysis and 

35 the change in the spectrum, which is very much different to 
the transient analysis of the present invention which is 
based on a direct transient detection in the time domain. 
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SUMMARY OF THE INVENTION 

The present invention provides an approach which is different 
in principle from all known methods for processing signals, 
5 According to the invention it has for example been found that 
the signal information relevant to recognition of speech is 
present in a transient part of the speech signal. Thus, the 
method of the present invention may involve a separation of 
the transient part of an auditory signal, a generation of a 
10 transient pulse corresponding to the transient part, and 
analysis of the shape of the pulse. In an auditory signal, 
the corresponding transient pulse may be repeated with time 
intervals, and the time interval of these periodic transient 
pulses is normally also analysed or determined. 

15 

In real life, the human ear reacts to energy changes at high 
frequencies in order to recognise phonemes or sound pictures. 
But in the present method transient pulses corresponding to 
the energy changes observed by the ear are extracted at these 

20 high frequencies, whereafter the transient pulses preferably 
are transformed to the low frequency range still maintaining 
the distinct features of the sound pictures or phonemes. 
Thus, by using the principles of the invention, it is 
possible to obtain distinct features within auditory signals 

25 by examining the transformed low frequency signals. 

The invention relates to the use of the shape of energy 
changes of a signal for identifying or representing features 
of the system generating the signal for example in 
30 recognition of sound features which can be perceived by an 
animal ear such as a human ear as representing a distinct 
sound picture are determined. 

The method of the present invention provides an expression 
35 for the transient conditions of the auditory signal. The 

method comprises a band-pass filtration of an auditory signal 
within the frequency range of the human ear and a detection 
of a low-pass filtered envelope, which envelope then can be 
analysed with known methods of signal analysis. The envelope 



2031 Fj nK 1 . go 1 / G H J / 3TM/ 13 l c i i i 9 




5 



is an expression of the transient part of the signal. 

The method of signal analysis, which should be used when 
analysing the envelope, and the characteristics of the band- 
5 pass filter, which should be selected, will depend on the 
purpose of the analysis. The purpose may be speech 
recognition, quality-measurement of audio products or 
acoustic conditions, and narrow band telecommunication. 

10 The invention also relates to a system for processing a 

signal to reduce the bandwidth of the signal with substantial 
retention of the information of the signal. The system may 
further comprise means for extracting the transient component 
of the auditory signal, and it may comprise means for 

15 detecting an envelope of the transient component. 

A signal may be separated into a sum of impulse responses 
generated by poles and zeroes in the system that has 
generated the signal, if the time between the excitation 
20 pulses are sufficient long compared to the duration of the 
impulse responses for the system. 

In WO 94/25958 it is shown that the envelope of the transient 
component in a speech signal is very important for its 
25 recognition and it is shown that the envelope of the impulse 
response will contain exponential functions and difference 
frequencies defined by the impulse response. 

A method based on damped sinus functions to extract important 
3 0 features from the envelope signal is described, and examples 
where the method is used on speech signals shows that the 
features are important in speech analysis . 

Before entering into a more detailed explanation of features 
35 of the method of the invention, a few definitions will be 
given; 

In short time analysis the transient component in a signal is 
a matter of definition. For auditory signals, the idea is to 
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obtain an expression that gives a response corresponding to 
the response in the cochlea to an abrupt change in the signal 
energy. An abrupt change in the signal energy corresponds to 
the transient component in the auditory signal. Thus, in the 
5 present context, the term "transient component" designates 
any signal corresponding to an abrupt energy change in an 
auditory signal. The transient component holds the signal 
information to be analysed and in order to analyse this 
information the transient component may be transformed to a 

10 corresponding transient pulse having a distinct shape. Thus, 
in the present context, the term "transient pulse" refers to 
a pulse having a distinct shape and substantially holding the 
information of the transient component of the auditory signal 
and thus corresponding to an abrupt change in the energy of 

15 the auditory signal. As mentioned above the transient part of 
a sound signal may be repeated with time intervals and thus, 
in the present context, the term "periodic" when used in 
combination with a transient component, response or pulse 
designates any transient component, response or pulse being 

20 repeated with intervals. 

The term "shape" designates any arbitrary time-varying 
function (which is time-limited or not time-limited) and 
which, within a given time interval Tp has a distinctly 

25 different amplitude level in comparison with the amplitude 
level outside the interval. Thus, Tp is the duration of the 
shape function when the shape function is time-limited, or 
the duration of the part of the function which has a 
distinctly different amplitude level in comparison with the 

30 amplitude level outside the time interval. 

In order to extract information from the shape of the energy 
changes, one broad aspect of the invention relates to 
represent the shape of the energy changes by the short time 
35 Laplace transform of a transient pulse of the signal. 

However, several methods can be applied in order to obtain a 
transient pulse corresponding to the change in energy, but it 
is preferred that an envelope detection is being used, where 
. the envelope preferably should be detected from a transient 
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response of the energy change in the auditory signal . 

The energy change representing the distinct sound picture can 
be a phoneme or vowel or any other sound which gives a sudden 
5 energy change in an auditory signal. 

It is also an aspect of the invention to provide a method for 
identifying, in an auditory signal, energy changes which can 
be perceived by an animal ear such as a human ear as 

10 representing a distinct sound picture, the method comprising 
comparing the shape of energy changes of the signal with 
predetermined energy change shapes representing distinct 
sound pictures. For the identification it is preferred that 
the shape of the energy changes are represented by the shape 

15 of a transient pulse of the signal, and it is furthermore 
preferred that the shape of the transient pulse should be 
obtained by an envelope detection of a transient response of 
the energy change in the auditory signal. 

20 The invention also relates to a method for processing a 
signal so as to reduce the bandwidth of the signal with 
substantial retention of the information of the signal, 
comprising extracting a transient part of the signal. The 
method may further comprise detecting an envelope of the 

25 transient part of the signal. 

Known methods of processing signals are based on a short time 
Fourier transform of signals, and it is assumed that the 
signals are steady state signals. 

30 

In steady state analysis the signal is assumed stable in the 
period the signal is analysed, and the steady state spectrum 
is calculated. 

35 In WO 94/25958 it is disclosed that transient pulses are 
important for speech coding and decoding in narrow band 
communication, for speech recognition and synthesis, and for 
sound quality in auditory products (i.e. loudspeakers, 
amplifiers and hearing aids) . 
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An important part of a transient signal is the exponential 
functions or damping ratios or time constants. The damping 
ratio is the reason that the impulse response has a finite 
5 duration. The fact that the transient signal is important for 
auditory perception indicates that the response from the hair 
cells are dependent of the time constants. If this is the 
case, it is possible that the damping ratios in the response 
from nerve cells in general are important for the human nerve 
10 system. 

Transient signals are also important in many other 
applications, among others signals generated by impacts from 
defects in rolling bearings and gear boxes. 

15 

Based on the transient signal, it is possible to determine 
the natural time constants and frequencies in the system 
generating the signal. Further it is possible to determine 
the excitation pulses of the system. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 



Fig. 1 shows a time-domain representation of a linear 
time- invariant system, 

25 

Fig. 2 shows the impulse response of a Butterworth low- 
pass filter of 3. order and a cut-off frequency at 
700 Hz, 

30 Fig. 3 shows the response with the filter relaxed for t< 
0 and with a 4000 Hz tone as input at t>0, 

Fig. 4 shows the s-plane with poles and the zero for 

35 

Fig. 5 shows H{a,co) for co^ and li;-, analysed parallel with 

the a axis, 



Fig . 6 
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shows transient characteristics in speech signals, 
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Figs . 



show processed speech signals, and 
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5 Fig . 13 shows a schematic of a filter bank according to the 



DETAILED DESCRIPTION OF THE DRAWING 

10 The importance of the transient part of a signal has been an 
overlooked phenomenon in signal analysis . 

The response of a linear system to either an impulse or a 
step function is defined by its transient response 
15 properties . 

The relationship between the input and the output for the 
linear time-invariant system shown in Fig. 1 can be written 
as the convolution of the input signal and the impulse 
20 response of the system: 



If the system is initially relaxed, and the input signal Vi(/) 
25 is zero for ^ < 0 then the lower integration limit of Eq . ( 1 ) 
can be replaced with zero. Eq. (1) then shows the important 
role played by the impulse response in terms of the actual 
signal processing that is performed by the system. It states 
that the input signal is weighted or multiplied by the 
30 impulse response at every instant in time and, at any 
specific point in time, the output is the summation or 
integral of all past weighted inputs. 

The impulse response of a real system has a finite duration 
3 5 and the transient response has the same duration . Fig . 2 
shows the impulse response of a Butterworth low-pass filter 
of 3. order and a cut-off frequency at 700 Hz. Fig. 3 shows. 
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10 

the response with the filter relaxed for t< 0 and with a 
4000 Hz tone as input at />0. 

In many processes v,(0 will be a pulse with a short duration 
5 and v/(/) 0 before the next pulse will be generated. 

The Laplace transform of a signal v{t) is defined by 

CO 

L{s) = \v{t)Q-"dt (2) 

0 

10 

00 

0 

If v(r) is the impulse response h{t) for a system with 2 
complex poles 

15 

h{t) =g-(^'^>^o)'_{.g-(*^o->^^o)' ^ ^ > 0 { 3 ) 

and 0 for t <Q and s ^ -{aQ± jco^) . 
20 the Laplace transform is 



5-hao 



or 

25 



a +0-0+ jo) 

H{a, CO) = ^^^^ ^€0,)){a ^a^-^j{(D -cy o)) 



From Eq.(4) it is seen that for {a^co) {-'a^,±co^) , //(cr, fi)) ±oo . 

30 

This. is a well known phenomenon and a logical . consequence of 
this is as follows: 
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If the signal analysed is dominated by the impulse response 
of the system generating the signal, it is possible to 
determine the natural time constants and frequencies for the 
5 system. 

Fig. 5 shows a plot of H(a/Co) for co-co^ and co-co-,. 

Analysing a signal along or parallel with the jco axis will 
10 give a frequency profile for a given a. 

Analysing a signal along or parallel with the a axis will 
give a time constant profile for a given jco. 

15 If a signal has a time constant profile with significant 

variations for specific frequencies, the signal is transient 
dominated. Opposite if the signal does not vary significantly 
for any frequency, the signal is steady state dominated. 

20 A short time Laplace transform is defined by: 



in which vi is the signal, L is the transformed signal, a is 
25 a time constant, and o) is an angular frequency. 

It is not possible to calculate the short time Laplace 
transform in the same way as DFT in the discrete time domain 

because e"" and , a and b reel, are not orthogonal 
30 functions for a^h . 

The short time Fourier analysis in the analogue time domain 
is based on a filter bank method. In this paper an equivalent 
method will be developed for the Laplace transform. 

35 

From Eq. (1) and Eq. (3) : 
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0 

t 

+ \vi{t-X)e-^"-'''^'dX (6) 

0 

5 

where u\t) is the complex conjugate of 7j{i) and we have 

10 Re[L(a,a;,0] = K(0 

From Eq. (6) and Eq, (7) it is seen that filtering the signal 
v.(r) by a filter with the impulse response h{a,Q),t) with 2 
complex poles will represent the reel part of the short time 
15 L{a,coj) transform. 



If we let v.(/) be equal to the impulse response of a single 
pole we have 

20 

t 



0 

25 

{a-a^)^ j{co-co^) 



and from Eq. (7) we have 

30 
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or 

5 



^ ^ 2^(a-ao)(g-" cos(tur)-g-°»' cosjco.t)) 



+ ;^ ^3 — z (9a) 



v.(0 ^ g"^^((Q'-C7,)cos(6;oO-(a^-6Jo)sin(6j,0) 
2A: (cr-cJo)' +(fy-£Uo)' 



'^((g- ao)c os(ay) - (a; - fi;^ ) sin(tyr)) 

: 7: — Z (9b) 



Eq. (9) is not defined for (a, a;) = (0X0,^^0 ) but from (8) we have 
in this case 



w(0 = >fce-^"<'^-'"°^'j^/A 
0 

15 

and 

2 0 (0 = Ikte-''^' cos(6)j) (10) 



10 



and we have v^(/)->0 for / -> 00 . 

25 Eq. (9) shows that the gain is inversely related to ct-cTq and 
co-o)q, and when (<Jq,o)q) is far from (cr,co) and e'^-e"''*'' is 
small, v^(0«0. For (cro,a)o) <- (a,iy) v^(0 will have Eq.(lO) as 
the limit. It is not immediately to see if Eq, (9) has the 
maximum energy for (cTq , ty^ ) <- (a, a;) . 
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In the DC domain Eq. (9) can be written as 



v„(0 = 2A:i— ^ (11) 

5 

The maximum for v^(l) can be found as follows 



10 when 



and Eq. (11) will have the maximum for this value. 
15 It can be shown that t„^^-~ when o^cf^. 

When cr^cTo we will have the approximated maximum with ^ = 



20 v,(^) = 2^^^ — ^ (13) 



From Eq. (13) it can be shown that 



2 5 for a->cro 



In Eq. (11) represent the signal to be analysed and e'"^ 

the filter. Table 1 shows the result with a filter having a 
100 s"' and the signal varying from 1 to 10000 s"^ 

30 
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15 

It is not surprising that the convolution acts as a low-pass 
filter. The important fact is that the exponential function 
in the DC domain in some way acts as frequencies do in the 
frequency domain. 

5 

In table 1 v^, (r^J is the result of a convolution where the 
signal is differentiated. The result is, as expected, a high- 
pass filter 

10 If we look on Eq.(9a) without exponential functions it can be 
written as 



Vo(/)= (14) 

15 

it is seen that for a; — » oc- we will have -> 0 . 



o- : 100 s'' 




L 






1 

10 
100 
1000 
10000 


s 

0, 046516871 
0, 025584279 
0, 010000000 
0, 002558428 
0, 000465169 


0, 954548457 
0, 774263683 
0, 367879441 
0, 077426368 
0, 009545485 


0, 009545485 
0, 077426368 
0, 367879441 
0, 774263683 
0, 954548457 



20 Table 1 v^(/^)is given by Eq.{ll, 12) and 

normalised by a and 2k. v^,(/^^)is a 
convolution where the signal is 
differentiated and normalised by 
2k. 
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For co«o)q we will have 



2k(s\n(cof) - s\n(coj)) 



(15) 



5 



It can be shown that for co cOq we will have 



(16) 



10 



This result is as expected unstable. 

In transient analysis only the beginning of the' signal is of 
interest, and if cOo»i Eq.(14) will act as a band-pass 
15 filter. 

Speech processing is based on fast energy pulse generated by 
the vocal cords or by friction in the articulation channel 
weighted by the impulse response in the articulation channel. 
20 The rise time for the excitation pulses has to be sufficient 
faster than the rise time of the energy of the impulse 
response . 

The shape of energy pulses are important features in speech. 
25 If the time between the pulses is periodical it is voiced 
speech, and if not it is unvoiced speech. For some phonemes 
abrupt changes in the energy pulses are important. 

From WO 94/25958 it is known that the shape of the energy 
30 pulses are important for speech recognition, especially the 
leading edge. In the following a method to extract features 
will be developed based on an envelope detection. 

The convolution expressed in Eq. (9) can be regarded as a 
35 response from 2 poles in the articulation channel excited by 
an impulse. If cf^^ ^ a we have from Eq(9a) 
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^oiO = . , , (sin(6<) - sin(ft>oOj (17) 



5 

The envelope is defined as 



eit) = ^u\t) + u\t) 



10 where 

2/(0 = «/(/)* ^ 

is the Hilbert Transform. 
15 The envelope of Eq. (17) is then 



e'^ I 

^o(0 = I r V(sin(tyr) - s\n{o)Qt)y + {-cQs{cot) + cos(<a;oO)^ 



2 0 =1 r^2(l-cos(ty-cyo)0 



i r(l-4-cos((w-cyo)0) (18) 



25 



The approximation is legal because |cos((£U - cOq )0| ^ 1 



As expected the envelope has a component with the difference 
frequency of the 2 frequencies. 

30 



20;il5DKl.i;)01/OHJ/:^TH/l j y.hccv. ".-^Si 



18 



The conclusion is that we can expect to find damped 
difference frequencies in the envelope of the transient 
component . 

5 To detect the damped difference frequencies a filter bank is 
used. The features might be detected as a convolution between 
the transient pulse and the impulse response of the filters. 

In general form the impulse response can be written as 

10 

h{i) = ke-'' s\n{f{X)t ^ ^) 

Where a = /I and o) - f {A) . 

15 In the following analysis /(/l) = L5A, /: = = 1.5/1 , and ^=0 are 
selected and we have 



h{i) = \5Xe"' sin(L5A0 (19) 

20 

By selecting = 1.5a Eq.(19) will act as a band-pass filter 
with a low Q in relation to the frequencies. Other ratios co/o 
than 1.5 may be selected and it is presently preferred that 
the ratio (co/o) ranges from 0.5 to 2.5. The exponential 
25 function gives the advance that it acts like natural time 
window that ensure that the signal is natural damped. The 
value of the parameters are selected by studying rise times 
in important transient pulses and by experiments. 

30 Fig. 6 shows transient characteristics in speech signals. The 
top figure shows 50 ms of an ''a" in ''hard key" pronounced by 
a female. 

The second signal is a band-pass filtration of the speech 
35 signal. The band-pass filter is a Butterworth filter with 6 
poles and a band width from 215.0 to 3550 Hz. This frequency 
band contains important transient pulses in the sensitive 
frequency interval of the ear. 
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The third signal is a energy detection of the transient 
characteristics of the band-pass filtered speech signal. The 
detection is an envelope detection performed by means of a 
5 rectification and a low-pass filtration of the signal. The 
filter is a Butterworth filter with 3 poles and a cut-off 
frequency at 700 Hz. 

In WO 97/09712 a method for automatically detecting the 
10 leading edges is disclosed. The method uses the maximum slope 
of the leading edge as reference, and the point before the 
maximum slope where the slope is less than a given threshold 
(10-20 % of the maximum slope) the leading edge is defined to 
begin. 

15 

The transient (envelope) signal in Fig. (6) has a DC 
component, which does not contain any information. Therefore 
it is preferred that the signal is differentiated before it 
is analysed e.g. by the filter bank shown in Fig. 13. 

20 

In Fig. 13, the filters (hi (t) , h2(t),..., hn(t)) in the filter 
bank connected between the input and the envelope detectors 
are band-pass filters having bandwidths corresponding to the 
bandwidths of the band-pass filters of the cochlea and having 
25 centre frequencies ranging from 1400 Hz to 6500 Hz. 

The output signals Oij (p) from the filter bank shown in Fig. 
13 is calculated by: 



30 



i=0, 1,..., 



N-1 



j = 0, 1,,.., 



M-1 



p < 0 



35 



p=0, 1,..., P-1 
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m=0, 1, ...,M-1 and M is the number of band-pass filters with a 
low Q in the filter bank connected between the outputs and 
the envelope detectors, p = 0,1,..., P-1 is the sample number, 
5 t' is the differentiated transient signal, and is the 

filter bank parameter and it is normalised by the sampling 
frequency. 

In the analysis M is selected to 10 and 1500 < < 12000 s~\ 
10 is not normalised. By this we have 1885 < OJ^ < 18850 s"^ or 

300 < f < 3000 Hz - 

This filtering process is not done in the cochlea but in the 
hair cells or • in the nerve system behind the hair cells. 

15 

The Figs. 7, 8,9, 10, 11, and 12 show the output of the 
processing of transient signals in the vowels, '^a'', ''o", '"i" 
in '"hard key" and ""soft key'' pronounced by a female and a 
male. Further the figures show plots of maxima of the output 
20 signals as a function of the time constant a of the 
corresponding filter . 

The figures show that maximum curves are very much alike for 
the same vowels, independent of whether it is pronounced by a 
25 female or male. 

With a library of templates and a distance measure it is 
possible to identify the sound picture, and it can be used 
for speech recognition and narrow band communication. 

30 

Thus, according to the invention a method and an apparatus 
are provided for determination of a parameter of a system 
generating a signal containing information about the 
parameter, in which the signal is short time transformed 
35 substantially in accordance with 




in which Vi is the signal, L is the transformed signal, a is 
a time constant, co is an angular frequency, and cp is a phase, 
or, in accordance with another transformation which will give 
5 rise to an L' (a,co,t) which in time intervals within which 
L(a,co,t) is larger than 10% of its maximum value is not more 
than 50% different from the result given by the short time 
Laplace transformation . 

10 In narrow band communication the transient pulses have to be 
identified and coded, and the decoder will contain a library 
of filters with corresponding transient response. The decoder 
library could also the transient responses. 

15 The present invention also relates to measurement of 

mechanical vibrations e.g. when testing devices that generate 
mechanical energy during operation, such as mechanical 
devices with moving parts, such as compressors for 
refrigerators, electric motors, household machines, electric 

20 razors, combustion engines, etc, etc. 

For example, it is known that measurement of vibration gene- 
rated or sound emitted by a device during operation can be 
useful for detection of malfunction of the device . Certain 
25 failures may generate sound or vibration of specific 
characteristics that can be recognised. 

The method may also comprise steps of classification for 
classifying a tested device in accordance with the determined 

30 parameters into one class of a set of predefined classes. 
Each predefined class may be defined by a* set of upper and 
lower limits for specific parameters determined according to 
the method. A device may then be classified as belonging to a 
certain class if its corresponding parameter values lie 

35 within corresponding upper and lower limits of the class. 
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Each class may correspond to a specific type of failure of 
the device. For example, shaft imbalance, wheel imbalance, 
crookedness, imperfections of teeth in cogs, tight bearing, 
5 loose bearings, etc, may cause the device to vibrate in 
different characteristic ways, whereby a characteristic 
mechanical vibration or sound is generated for each type of 
failure. The type of failure of the device may then be 
detected by comparing determined device parameters with 
10 corresponding parameter values of various predetermined 
classes. 

The upper and lower limits of a specific class of devices may 
be determined by testing a set of devices known to belong to 
15 that class. For example, the upper limits may be determined 
as the average of specific parameter values plus three times 
the standard deviation. Likewise, the lower limits may be 
determined as the average of parameter values minus three 
times the standard deviation. 



20 
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CLAIMS 

1 . A method for determination of a parameter of a system 
generating a signal containing information about the 

5 parameter, comprising the step of short time transforming the 
signal substantially in accordance with 

0 

in which Vi is the signal, L is the transformed signal, a is 
10 a time constant, co is an angular frequency, and cp is a phase. 

2. A method according to claim 1, wherein the step of 
transforming comprises filtering the signal Vi with a filter 
having a pole at a + jcot and a pole at a - jcot. 

3. A method according to claim 1 or 2, comprising steps of 
transforming the signal Vi for a plurality of sets of a and co 
values • 

20 4. A method according to any of the preceding claims, further 
comprising the step of determining a maximum of at least one 
transformed signal L(a,a),t). 

5. A method according to any of the preceding claims, further 
25 comprising the step of comparing transformed signals L with 

corresponding reference signals in order to determine 
parameters of the system. 

6. A method according to any of the preceding claims, further 
30 comprising a step of pre-processing the signal before the 

step of short time transforming, the pre-processing being 
selected from the group consisting of filtering, 
rectification, differentiation, integration, and 
amplification . 

35 
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7. A method of transmitting a signal containing information 
of a set of parameters of a system generating the signal, 
comprising processing the signal according to any of the 
preceding claims and further comprising the step of 

5 transmitting the determined parameter values, 

8. A method according to claim 7 further comprising the step 
of regenerating the signal from the determined and 
transmitted parameter values. 

10 

9. A method of transmitting a signal containing information 
of a set of parameters of a system generating the signal, 
comprising processing the signal according to any of the 
preceding claims and further comprising the steps of 

15 

comparing the signal with a library of signals generated for 
a predetermined set of parameter values by the system, 

selecting the library function that constitutes the best 
20 match to the signal, and 

transmitting an identification signal that identifies the 
matching library function. 

25 10- A method according to claim 9, further comprising the 
steps of receiving the identification signal and generating 
the corresponding library signal. 

11. A method of classifying a system according to one or more 
30 parameters of the system generating a signal containing 
information about the one or more parameters, comprising 
determining the one or more parameters according to any of 
claims 1-6 and further comprising the step of classifying the 
system in accordance with the one or more determined 
35 parameters into one class of a set of predefined classes 

defined by predetermined ranges of values of the parameters. 
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12. A method for communicating an auditory signal, comprising 
processing the signal by the method according to any of 
claims 1-6, transmitting the processed signal, and receiving 
the processed signal by a receiver. 

5 

13. A method according to claim 12, wherein, prior to 
transmission of the processed signal, the signal is coded 
into a digital representation, and the coded signal is 
decoded in the receiver so as to reestablish transient pulse 

10 shapes perceived by the animal ear such as the human ear as 
representing the distinct sound pictures of the auditory 
signal . 

14 . A method according to claim 13, wherein the digital 

15 transmission is performed at a bandwidth of at the most 4000 
bits per second. 

15. A method according to claim 14, wherein the bandwidth is 
at the most 2000 bits per second . 

20 

16. A method according to claim 15, wherein the bandwidth is 
in the interval of 800-2000 bits per second. 

17. A method according to any of claims 13-16, wherein a 

25 second and further pulses in a sequence of identical pulses 
are represented by a digital value indicating repetition. 

18. A method according to any of claims 1-6, comprising 
filtering the signal Vi in a filter bank comprising a 

30 plurality of band-pass filters interconnected in parallel 

with centre frequencies ranging from 1400 Hz to 6500 Hz, each 
of which is connected in series with an envelope detector and 
a filter bank comprising a plurality of low-pass filters 
interconnected in parallel and having cut-off frequencies 

35 ranging from 300 Hz to 3000 Hz and time constants a ranging 
from 1500 s"^ to 12000 s~\ 
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19. An apparatus for determination of a parameter of a system 
generating a signal containing information about the 
parameter, comprising "a processor that is adapted to short 
time transform the signal substantially in accordance with 



in which Vi is the signal, L is the transformed signal, o is 
a time constant, co is an angular frequency, and cp is a phase. 

10 20. An apparatus according to claim 19, wherein the processor 
comprises a filter for filtering the signal vi and having a 
pole at a + jcot and a pole at o - jot. 

21. An apparatus according to claim 19 or 20, wherein the 
15 processor comprises a plurality of filters for filtering the 
signal Vi, each filter having a different set of a and o) 
values. 
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3. Order. LP, 700 Hz, Butterworth 



impulse response. 
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3. Order. LP. 700 Hz. Buttervvorth 
Step frequency response. 
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50 ms 




Speech signal 




Transient isolation in the speech signal, band width 2150-3550 Hz 




Energy detection of the transient pulses by means of envelop detection, 
rectified and low pass filtered at 700 Hz 
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